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Abstract 

Di-jet event rates have been measured for deep-inelastic scattering in the kinematic 
domain 5 < < 100 GeV^ and 10~^ < xbj ^ 10~^, and for jet transverse momenta 
squared Pt ^ Q^- The analysis is based on data collected with the HI detector at 
HERA in 1994 corresponding to an integrated luminosity of about 2 pb~^. Jets are 
defined using a cone algorithm in the photon-proton centre of mass system requiring 
jet transverse momenta of at least 5 GeV. The di-jet event rates are shown as a 
function of and xbj- Leading order models of point-like interacting photons 
fail to describe the data. Models which add resolved interacting photons or which 
implement the colour dipole model give a good description of the di-jet event rate. 
This is also the case for next-to-leading order calculations including contributions 
from direct and resolved photons. 
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1 Introduction 



The study of jets in deep-inelastic lepton-proton scattering (DIS) provides a testing ground 
for perturbative QCD. Partons emerging from the scattering process manifest themselves 
as jets of coUimated hadrons which are observable in the experiment. 

In the naive quark parton model, the virtual photon is absorbed by a single quark 
(antiquark) of the proton resulting in one jet from the struck quark and one from the 
proton remnant. Both jets have no transverse momentum in the photon-proton centre 
of mass frame (cms), when neglecting the intrinsic motion of the partons inside the pro- 
ton. To first order in a^, the leading order (LO) for di-jet production, two jets with 
balanced transverse momenta in the photon-proton cms are produced in the hard scat- 
tering process, in addition to the proton remnant jet. The hard scattering can either be 
the quark initiated QCD-Compton (QCDC) or the gluon initiated photon-gluon fusion 
(BGF) process. 

In this analysis we present a measurement of the fraction R2 of di-jet events in all 
DIS events, referred to as the di-jet rate. It is presented as a function of the Bjorken 
scaling variable xbj, integrated over the virtuality of the exchanged photon Q^, and of 
Q^, integrated over xbj- Jets are defined using a cone algorithm in the photon-proton cms 
requiring jet transverse momenta of at least 5 GeV. The measured jet rates are corrected 
for detector effects. 

Previous measurements of jet rates at HERA |^ used the JADE jet algorithm at 
photon virtualities large compared to the squared jet transverse momenta pf. The 
present analysis probes a region of jet phase space characterized by jet transverse momenta 
squared of similar size or larger than the photon virtuality, Pt/Q"^ ^ ^- It has significantly 
better precision and has its emphasis on higher compared to a previous study of single 
inclusive jet production [0]. There the data were found to be in good agreement with 
LO QCD models which included a resolved partonic structure of the virtual photon that 
evolves with Q^. 

In this study we investigate whether the di-jet rate can be described by LO QCD 
models with just point-like (direct) interactions of the virtual photon and with models 
with additional contributions from resolved photons, which may be considered as an 
effective description of higher order QCD effects. We also consider the colour dipole 
model. Finally our measurements are compared to next-to-leading order (NLO) QCD 
calculations which include either only direct or direct and resolved virtual photons. 

The measurement was performed using data taken in 1994 with the HI detector at 
the HERA storage ring, where 27.5 GeV positrons were collided with 820 GeV protons. 

2 The HI Detector 

A detailed description of the HI apparatus is given elsewhere 0. The parts of the detector 
which are essential for this measurement are the liquid argon (LAr) calorimeter Q], the 
backward lead-scintillator calorimeter (BEMC) and the tracking chamber system. 
The energy of the scattered positron is measured in the BEMC which covers the range 
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in polar angleQ 9, from 151° to 176°. It consists of stacks of lead and scintillator plates with 
a total of 21.7 radiation lengths. The BEMC is laterally segmented into square modules 
of 16 X 16 cm^, with smaller modules at the inner and outer radii. The scintillation light 
is read out with photodiodes via wave length shifters along two opposite sides of each 
module. The absolute energy scale was determined to a precision of 1% The energy 
resolution is given by ctb/E = 39%/E © 10%//E © 1.7% {E in GeV) §. 

A cluster energy deposition exceeding a threshold of ~ 7 GeV in the BEMC was the 
primary trigger condition for events used in this analysis. 

The position of the scattered positron is measured with the backward proportional 
chamber (BPC) located in front of the BEMC covering the angular range 155° < 9 < 
174.5°. The BPC consists of four layers of wires strung vertically, horizontally, and at 
±45°. The position resolution is ax,y = 1.5 mm. 

Hadronic energy is detected in the highly segmented (~ 45000 channels) LAr calorime- 
ter which extends from 4° < 9 < 154°. The depth of the LAr calorimeter varies between 
4.5 and 8 hadronic interaction lengths in the region 4° < ^ < 128°. The uncertainty of the 
energy scale for hadrons is 4%. The hadronic energy resolution is as/E = 50%/^/E(B2% 
{E in GeV), as measured with test beams [0. 

Charged tracks in the central region (25° < 9 < 155°) are measured with the central 
drift chamber system. Two jet chambers with wires in the z-direction allow measurements 
of track positions in the r-0-plane to a precision of cxr^ = 170 /im. The z coordinate is mea- 
sured to a precision of o"^ = 320 /xm using drift chambers with wires forming approximate 
circles around the beam. The momentum resolution is (Jpjp^ < 1% GeV~^. 

The forward tracking detector covers 7° < 9 < 25° and consists of drift chambers with 
alternating planes of parallel wires and others with wires in the radial direction. It allows 
the measurement of track segments to a precision of a^^y < 200 fim. 

Two electromagnetic calorimeters located downstream in the positron beam direction 
measure positrons and photons from the bremsstrahlung process ep epj for the purpose 
of luminosity determination. 

3 Data Selection 

The data sample used for the present analysis corresponds to an integrated luminosity 
of 1.97 pb~^ taken during the 1994 run period. The phase space region of DIS events 
considered in this analysis is defined as follows: 

156° < 9e < 173° 

E' > 11 GeV (1) 
y > 0.05 

Here, 9^ is the polar angle of the scattered positron and E' is its energy. The variables xbj, 
and y (the inelasticity variable) are all determined from the 4-vector of the scattered 
positron. This selection ensures that the scattered positron is well inside the acceptance 
region of the BEMC, that the trigger efficiency is high, that the kinematic variables 

^The polar angle 9 is defined with respect to the positive z-axis, the proton beam direction. 
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are well reconstructed, and that photoproduction background and radiative corrections 
are small. Photoproduction events, where the scattered positron is not detected in the 
backward direction, form a background if a particle from the hadronic final state entering 
the BEMC is misidentified as the scattered positron. 

Additional cuts are applied for the identification of the scattered positron and to 
further suppress the influence of QED radiation and photoproduction background |^: 

• The event must have a reconstructed vertex with a z position within ±30 cm of the 
nominal position. 

• The candidate positron shower is required to have a small lateral spread by ap- 
plying the cut Tciust < 5 cm, where r^ust is the energy- weighted mean transverse 
distance from the shower centre of gravity of each energy deposition sampled by the 
photodiodes. 

• There must be a BPC signal within 5 cm of the straight line connecting the shower 
centre with the event vertex. 

• The quantity J2i{Ei —pz,i), where the sum is over all calorimeter energy depositions 
in the final state, is expected to be equal to twice the positron beam energy. An 
undetected positron in a photoproduction event or initial state photon radiation will 
decrease the value of this observable. For this analysis, 35 < J2i{Ei —pz,i) < 70 GeV 
is required. 

4 Jet Reconstruction and Selection 

Jets are reconstructed using clusters of energy measured in the LAr calorimeter. Clus- 
ter energies are corrected for the difference in response to hadronic and electromagnetic 
energy deposition and for losses due to dead material and cracks. The cluster energy 
and the direction from the interaction point to the cluster centre are used to construct a 
massless four-vector. 

The calorimetric energy measurement can be improved for low energy particles by 
using in addition to the energy the measured momentum of each charged particle track. 
To avoid double counting of energy each track was allowed to contribute at most 300 MeV. 
The value of 300 MeV was found to be optimal for reconstructing the transverse momenta 
of jets in simulated events [§. 

Jets are defined in this analysis using a cone algorithm |p. A cone is defined by a 
circular area of radius R in the r]* - 0* plane, where rj* and 0* are the pseudo-rapidity^ 
and azimuthal angle in the photon-proton cm^. A jet candidate consists of all objects 
(clusters and tracks) whose massless four vectors fall inside a cone. The jet transverse 
momentum Ptj^^t is calculated as the scalar sum of the transverse momenta of the jet 
objects. The jet rj* and (p* (jet direction) are calculated as the p^-weighted averages of 
the rj* and 0* of the objects. This way of calculating the jet parameters is usually called 

^The pseudo-rapidity rj* is given by — lntan(6'*/2). 

•^Variables measured in the photon-proton cms are marked by a * . The positive z direction is defined 
to be along the incident proton direction. 
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the "pi"-scheme [|TU|. An iterative procedure is used to find the jets of an event. Initially, 
every object in turn is used to define the cone centre of a candidate jet. The jet directions 
of the candidate jet^ are then used as the cone centres for the next iteration. This is 
repeated until the resulting jet directions are identical to the cone centres. Then, also the 
midpoint in the f]* - 0* plane of each pair of jets is considered as a candidate jet centre, 
and the procedure is repeated. Jets which have more than a fraction / of their Pt^^t 
contained in a higher transverse momentum jet are discarded. Finally, Pt^^t is required to 
exceed a minimum value Pf min- 

In this analysis, the following parameters were chosen: 

R=l, pI^.^ = 5 GeV, and / = 0.75. (2) 

Exactly two jets per event fulfilling these criteria are demanded. In addition, the pseudo- 
rapidity difference Arj* of the two jets is required to be in the range 

\Ar]*\ < 2. (3) 

In leading order this cut is equivalent to requiring | cos 6*1 < 0.76, where 9 is the polar 
angle between the emerging and incoming partons in the parton-parton or gamma-parton 
cms. It separates the jets from the proton remnant. The resolution in jet transverse 
momentum Ap^/p^ is approximately 20% at ~ 5 GeV. 

The number of di-jet events found is 4 957 while the total number of DIS events selected 
amounts to 112 806. To obtain the di-jet rate R2, the number of di-jet events is divided 
by the total number of events in the same region of XBj and Q^. R2 is measured in bins 
of XBj, integrated over Q"^, and in bins of Q^, integrated over xbj- 



5 Data Correction 

The residual background from photoproduction processes was determined using the PHO- 
JET Monte Carlo (MC) generator [|rT| and was separately subtracted from the total num- 
ber of events and the di-jet events as a function of and xbj- This generator has been 
proven to give a good description of photoproduction background The correction for 
this background as well as the other corrections described below were obtained using MC 
events which were processed by the HI detector simulation, reconstruction, and analysis 
chain. The largest subtraction of the photoproduction background occurs in the lowest 
XBj and bins, where it amounts to 14% and 9% of the total event sample respectively, 
and to 3% and 1% of the di-jet sample. It is below 5% in the total sample and negligible 
in the di-jet sample in all other bins. For the di-jet rate, the correction is only significant 
in the lowest Xsj and bins where it increases the rate by ~ 10%. 

Radiation of photons from the incoming or outgoing positron leads to values of XBj, Q^, 
and y, as determined from the scattered positron, which differ from the true kinematics 
of the photon-proton interaction.^ These effects are different for the total and the di-jet 
sample. They were corrected using the DJANGO MC generator |jl2|. The correction 
factor on the di-jet rate was found to be 1.08, independent of XBj and Q^. 

"^Several initial cone centres may result in the same candidate jet. 

^In the outgoing positron case, this only applies if the angle between the radiated photon and the 
outgoing positron is large. 
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Figure 1: Transverse energy flow with respect to the jet axis, a) versus f]* — t]*^^ in a 
slice defined by |0* — < 1, b) versus 0* — 0*,^ in a slice defined by \ri* — r]*,^ < 1. 
In the Figs, c) and d) the rj* and spectra of the jets are shown. In Fig. e) the \Ari*\ 
distribution is displayed. Points denote data, the histograms indicate the distributions 
obtained from LEPTO (dashed) and ARIADNE (dotted). The curves are normalized to 
the number N of entries; there are two entries per di-jet event in Figs. a)-d) and one in 
Fig. e). For Fig. d), Pt^min (Eqn. ^ was lowered to 3.5 GeV and for Fig. e), the cut on 
|A?7*| (Eqn. ^ was omitted. 
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The correction of the di-jet rate for detector acceptance and efficiencies was performed 
with two MC models, LEPTO ||13| and ARIADNE , which will be discussed in section 



7. They were used with two sets of parton density parameterizations, MRS-H ||T5[ and 



GRV-94 HO [0 as implemented in PDFLIB ||T^. The average prediction of these models 
was used to obtain bin-wise correction factors c = ]^^c,hadronYj^MC,recon. di-jet rate 

R2. This procedure is justified since the observables of the jet events, which are sensitive 
to detector effects, are well described by the Monte Carlo simulations. Critical observables 
in this sense are the energy fiow within and around the jets, the rj* and distribution of 
the two jets and the pseudo-rapidity difference |A?7*| between them. 

Fig. |l] shows a comparison of the experimental distributions with the two MC models. 
Only the curves obtained with the MRS-H parton density are shown, the curves with 
GRV-94 HO are very similar. There is good agreement between data and simulated 
events except for the 77* distribution of the jets. On top of the bin-wise correction, which 
ignores this small discrepancy, an additional correction for this effect was applied 0. It 
takes into account the fact that on average the 77* of the jets for data is higher than for 
MC which leads to an overestimation of the correction factors c, as the jet reconstruction 
efficiency depends on rj*. This correction reduces R2 by 10% in the lowest and by 1% in 
the uppermost bins of and xbj- The combined correction factors vary between 1.0 for 
low XBj and values and 1.2 for high values. 



6 Systematic Errors 

Several sources of systematic uncertainties were investigated. A change in the hadronic 
energy scale of the LAr calorimeter by its estimated precision of ±4% results in a global 
change of the di-jet rates by The correction for radiative effects has a global uncer- 
tainty of ±3% based on Monte Carlo statistics. These two errors were added in quadrature 
to give an overall systematic error of ^gy^° . 

Changing the energy scale of the positron measurement in the BEMC within its ±1% 
uncertainty results in a change of R2 by ±2% in all bins. In the xbj bins, the change 
varies between ±1% for the lowest and ±9% for the highest bin. The systematic errors on 
the corrections for acceptance and efficiency were obtained by using the maximal variation 
of the correction factor for any particular model compared to the mean in each bin. They 
are of the order of 10%. The additional corrections for the difference in the mean values of 
rj* between experiment and simulation have a systematic error of the order of 2%. These 
errors were added in quadrature to give a bin by bin systematic error. It varies between 
5% and 19%. 



7 QCD Calculation of Di-jet Rates 

Scattering processes involving the production of high partons (hard scattering pro- 
cesses) are expected to be well described by perturbative QCD. In this analysis hard 
collisions are selected by requiring two jets with transverse momentum above 5 GeV. We 
consider three different QCD inspired models labeled DIR (direct), DIR+RES (direct -|- 
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a) direct 



b) resolved 



c) CDM 




Figure 2: Generic diagrams of initial state parton emission in ep scattering (a,h). In the 
direct process (a) the hardest emission given by the QCD matrix element occurs at the top 
of the ladder. Emissions down the ladder are ordered with decreasing transverse momenta 
kt. In the resolved process (b) the hardest emission given by the QCD matrix element may 
occur anywhere in the ladder with increasingly soft emissions along the ladder towards the 
proton and the photon. In the colour dipole model (c) gluon emissions are not ordered in 
transverse momentum kt. 

resolved), CDM (colour dipole model), and NLO QCD calculations at the parton level for 
comparison with experimental data. We briefly describe their most important features. 

DIR: The LO QCD matrix elements (BGF and QCDC) are convoluted with the parton 
densities in the proton. Only direct interactions of the photon are considered as 
indicated in Fig. ^a. For this model we have used the LEPTO [|13| and RAPGAP 
[0 MC programs. The latter has been used without generating diffractive (rapidity 
gap) processes. RAPGAP and LEPTO give results consistent with each other to 
better than 10%. 

DIR+RES: In addition to the direct contribution discussed above, a contribution from 
quarks and gluons inside the photon is considered, as shown in Fig. ^p. This resolved 
photon contribution is assumed to set in only for scales of the hard subprocess ^Q^- 
For the virtual photon parton densities the SaS-lD parameterization was used. 
For this set of parton densities, the scale of the onset of the anomalous contribution 
in the virtual photon Pq = max(QQ, Q"^) was chosen, where Qq = 0.36 GeV^ is 
the starting point of the evolution]^. The LO resolved photon contribution is 
implemented in the RAPGAP program. The contribution from longitudinal virtual 
photons is neglected. 

^this corresponds to the SaSgam parameter IP2=2 
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In both the DIR and DIR+RES models we have used for the renormahzation and 
factorization scale /i^ = + pi as it provides a smooth transition between the DIS 
and photoproduction regimes. Additional emissions in the initial and final state are 
generated by parton showers in the leading log DGLAP |21| approximation. In this 



approximation the radiated partons in the initial state are strongly ordered in transverse 
momentum /c^Q, with the hardest emission in the ladder occurring next to the hard matrix 
element (Figs. |^ and ^jb). 



CDM: In the colour dipole model as implemented in the MC generator ARIADNE 
, gluon emission originates from a colour dipole stretched between the scattered 
quark and the proton remnant. Each emission of a gluon leads to two dipoles 
which may radiate further, generating a cascade of independently radiating dipoles 
(Fig. ^). These gluons are not ordered in kt. A similar feature is found in the 
BFKL evolution scheme. 

The colour charge of the proton remnant (a di-quark in the simplest case) is assumed 
not to be point-like, leading to a phenomenological suppression of gluon radiation 



IJ] in the direction of the remnant. This suppression occurs for hard gluons with 
wavelengths smaller than the size of the remnant. In addition, the colour charge of 
the scattered quark is taken to be extended, depending on the virtuality of the 
photon (photon size suppression). This in turn leads to a suppression of radiation 
in the direction of the scattered quark [|14|, p^ . 

The QCDC component of the di-jet rate depends in the CDM model on the size 
of the colour charge while for the DIR model it depends on the parton densities of 
the proton. At low xsj this results in a considerably enhanced di-jet rate for CDM 
compared to DIR [ESI. The photon-gluon fusion process, which is not naturally 



described by the CDM, is treated similarly as in the DIR approach discussed above. 

In the DIR, DIR+RES, and CDM models, hadronization was performed with the Lund 
string fragmentation scheme as implemented in JETSET [2^. 



NLO: Finally, we consider two calculations in next to leading order (NLO) in the strong 
coupling constant Us as implemented in the Monte Carlo integration programs DIS- 
ENT and JETVIP [|28|. These programs provide cross sections for partons 



rather than a full hadronic final state. DISENT takes the soft and coUinear diver- 
gencies arising in any NLO QCD calculation into account by using the subtraction 
method while JETVIP relies on the phase space slicing method. Both DISENT 
and JETVIP calculate NLO cross sections assuming a direct interacting photon. In 
addition, JETVIP provides a consistent calculation in NLO of direct and resolved 
interacting photons using parameterizations of the virtual photon structure func- 
tions. For the latter, the SaS-lD parameterization |19| (transformed to MS) was 
used. For both the DISENT and JETVIP calculations the factorization and renor- 
mahzation scales were chosen to be /x^ = + 50 GeV^, where 50 GeV^ represents 
a good estimate of the average transverse momentum squared of the jets in the 
hadronic cms for the selection described before. 



^kt is the transverse momentum relative to the proton and photon axis in the photon- proton cms. 
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For comparing the corrected di-jet rate with models and parton level calculations, we have 
used the CTEQ4M parameterization ||29| of parton densities inside the proton with the 
corresponding of 202 MeV (different parton density parameterizations were used in 
the models used to correct for detector effects, c.f. section 5). The DIR, RES, and CDM 
models implement the one loop expression for the calculation of and in the DISENT 
and JETVIP programs the two loop expression was used. 



8 Results and Discussion 

The di-jet rate R2 is shown in Figs. |^ and |b as a function of and x^j respectively. The 
data have been corrected for detector effects to the hadron level. The results correspond 
to the phase space region defined by Eqs. The data show a jet rate rising with 

and flat in xsj except for the highest x^j-value. 

For reasons to be explained later, two further scenarios have been investigated, where 
in addition to our basic requirement of pj' > 5 GeV for each jet (symmetric scenario), we 
demand either at least 7 GeV for the jet with the highest pi (asymmetric scenario) or 
at least 13 GeV for the sum of the absolute values of the jet transverse momenta (sum 
scenario). The results for R2 for the asymmetric scenario are given in Figs. ^,d and for 
the sum scenario in Figs. |^,f. Table 1 in the appendix summarizes the di-jet rates for 
the three different scenarios. 



8.1 Comparison of data with LO QCD models 

The results for the three different selections of jet phase space are compared to predictions 
from MC models based on perturbative QCD (see previous section). The LO DIR model 
fails to describe the data as demonstrated in Figs. ^ to ^f with RAPGAP. In particular 
in the region of small and xbj the DIR model underestimates the data by a factor 2-3. 

Choosing /x^ = for the hard scale does not change the results considerably. Using 
different parton density parameterizations (CTEQ4L, CTEQ4A4, CTEQ4HJ, MRSRl, 
and MRSR2) leads to variations in R2 of up to 10% in the lowest and the highest 
and XBj bins when compared to CTEQ4M, our default. The world average value of 



a,(M|) = 0.118 corresponds to AgL = 209li MeV 13, very close to the fit value of 



202 MeV for CTEQ4M. With the CTEQ4M parton densities but with Ag^ = 250 MeV 
(an increase of about one standard deviation), R2 increases by less than 10% in the lowest 
and Xsj bins. 

We conclude that a LO matrix element calculation assuming only direct interactions 
of the virtual photon in combination with DGLAP parton showers as an approximation 
to higher order effects is not able to account for the observed di-jet rates. 

Adding a significant contribution to the di-jet cross section from resolving the structure 
of the virtual photon, as predicted by the DIR-I-RES model as implemented in RAPGAP, 
gives a good description of the data. It should be noted, however, that considerable 
freedom exists in tuning the model to data, in particular by varying the choice of the 
hard scale, and the parton densities in the virtual photon. 
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Figure 3: Di-jet rate R2 as a function ofQ^ (a,c,e), integrated over xbj, and as a function 
ofxBj (b,d,f), integrated over Q"^ for the symmetric (a,b), asymmetric (c,d), and sum (e,f) 
cut scenario on the pi of the two jets. The data are corrected to the hadron level. The 
inner bars give the statistical errors, the full error bars include the bin by bin systematic 
errors. The grey band shows the overall systematic uncertainty due to the hadronic energy 
scale of the calorimeter and the uncertainty of the radiative correction. Also indicated are 
the predictions from MC models. Their statistical errors are smaller than the statistical 
errors of the data. 
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The CDM model, as implemented in ARIADNE, is also able to describe the di-jet rate 
well, both in absolute value and in the and xsj dependence. We used a parameter 
setting which had been tuned to give a good description of transverse energy flows and 
particle spectra ||31]]. Here too, it should be remarked that by varying the parameters 
for the proton and photon size suppression (see section 2) within sensible limits^ the 
predictions of this model can be changed by up to 40% in the lowest bin and about 20% 
in the highest bin in and Xsj- 



8.2 Comparison of data with NLO QCD calculations 

We now investigate whether a NLO QCD calculation is able to describe the data. For 
this purpose we have used results from the programs DISENT and JETVIP . For the 
calculation of R2 for the direct or point-like coupling of the photon to the partons in the 
proton, the two programs agree to better than 5%. As mentioned in section 7, JETVIP 
can also calculate the direct and resolved photon contribution in NLO. Both programs 
provide parton level cross sections rather than a full hadronic final state. However, the 
DIR and CDM models suggest that the hadronization effects are small for jet transverse 
momenta above 5 GeV. The di-jet rate at the parton level was found to be for LEPTO 
(ARIADNE) typically 9% (20%) and not more than 12% (25%) higher than the rate at 
the hadron level. 

Figs. 4a,b show the hadron level di-jet rate R2 versus and XBj for the asymmetric 
scenario, and in Figs. 4c, d for the sum scenario, compared to the NLO QCD calculation 
of the direct contribution (labeled DIR in the figures) by JETVIP and DISENT. Good 
agreement is observed between data and the direct NLO QCD calculation, except for the 
lowest and XBj bin. 

The sensitivity to variations of the parton density parameterizations is similar to the 
LO DIR case discussed above. Varying the factorization and renormalization scale /x^ by 
factors of 4 results in cross section variations of less than 20%. Choosing = as the 
scale enhances the cross section in the lowest xsj and bins by up to 30%, improving 
the agreement with data. At the same time this introduces however a large sensitivity 
to scale variations (up to 50% and 65% in the lowest and Xsj bins). This indicates, 
as one might expect, that is not the proper scale to use in a kinematic domain where 

« pI 

The agreement with data at low xsj and is improved when contributions from 
resolving the virtual photon structure are included in NLO (labeled DIR-PSP+RES in 
the figure). In order to avoid double counting in the full NLO QCD calculation it is 
necessary to subtract the contribution from the virtual photon splitting into gg, where one 
of the quarks subsequently interacts with a parton from the proton to produce two high pt 
jets, since this contribution is part of the parameterization of the virtual photon structure 
function ||T9l . We refer to this perturbatively calculated contribution from photon splitting 



as defined in [^] as PSP and the contribution from resolving the photon structure as RES. 

Two interesting observations can be made. First, at large the difference between 
the NLO direct part (DIR) and the full calculation (DIR-PSP+RES) is found to be rather 

^PARA(IO) and PARA(15), default 1.0 Q; we used 1.5 and 0.5 respectively and varied them inde- 
pendently between 0.5 and 1.5. 
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Figure 4: Di-jet rate R2 as a function ofQ"^ integrated overxsj, and as a function of 

XBj (b,d), integrated over Q"^ . The data (corrected to the hadron level) for the asymmetric 
(a,h) and the sum scenario (c,d) are compared to different NLO calculations (at the parton 
level). The data are the same as those shown in Figs. 3c to 3f. The statistical errors of 
the NLO calculations are smaller than the statistical errors of the data. 
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small, which implies that the NLO resolved part (RES) is saturated by the contribution 
from virtual photon splitting (PSP). Second, the full NLO calculation is close to the LO 
RAPGAP DIR+RES prediction shown as the full line in Figs. 3c to 3f. This suggests, 
together with the first observation, that the large resolved contribution needed in LO to 
describe the data for the larger bins is included in the NLO DIR cross section. 

It should be noted that in NLO the RES contribution depends less on the choice of the 
hard scale and the parton densities in the virtual photon than in LO. This is due to the 
subtraction procedure and because it is a NLO calculation . Of course the uncertainty 
due to the rather poorly known parton density of the virtual photon remains. 

A comparison of the data on R2 and the NLO QCD calculation for the symmetric 
scenario is not shown, because the calculation for this case is not reliable although the 
measurement is valid and infrared safe. The calculations from both DISENT and JETVIP 
underestimate the data and give different predictions. This can be understood as a feature 
of any fixed order calculation which gives large negative cross sections in the phase space 
region where both jets have almost identical transverse momenta. The problem in the 
prediction of di-jet rates for symmetric cuts on the jet transverse momenta has been noted 
in the framework of the phase space slicing method and discussed in detail for this 



and the subtraction method in ||33[ . A correct treatment of this phase space region would 



need a resummation to all orders 33, 34 



8.3 Event topology 

The conclusions of the underlying picture derived above can be checked by a study of 
the event topology. In the DIR+RES model the hardest emission leading to the observed 
jets may occur anywhere in the ladder, as depicted in Fig. ||b. In this case, additional 
hadronic activity is expected from the virtual photon "remnant" and additional parton 
emission from the top part of the ladder. 

This activity is expected in the direction of the virtual photon which corresponds to 
the backward region of the detector. Similar hadronic activity in the backward region 
can also be expected from the CDM model due to the absence of kt ordering between the 
photon and the proton vertex (see section 2). Both models predict a dependence of 
this effect which increases as approaches zero. 

We define an observable which is sensitive to additional energy flow in the photon 
direction: 

had. final state 

In the limit ^ and in LO this corresponds to the fractional momentum of the parton 
from the photon entering the hard subprocess and giving rise to the observed jet system. 
In this picture 1 — x'^'^ corresponds to the fractional energy of the photon remnant. 

Fig. I^a shows the uncorrected distribution of x'^^ for data in three different ranges 
of Q^. The MC events which were used for comparison in Fig. ^ have been subject to a 
detailed simulation of the HI detector. In the data an increase at low x'^^ is noticed as 
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Figure 5: Uncorrected distribution of x"^^ in three different bins (a). In b) to d) 
the data are compared for each bin in to the DIR model and the DIR+RES model as 
given by RAPGAP, and the CDM model as implemented by ARIADNE. The figures are 
normalized to the number N of di-jet events. The error bars indicate the statistical error 
only. 



18 



decreases. No such effect is seen for tlie DIR model as represented by RAPGAP (full 
line in Figs, to ^). The DIR+RES model of RAPGAP is able to give a reasonable 
description of both shape and dependence (dashed line in Figs, to ^). The CDM 
model shows a similar dependence but fails to describe the shape of the distribution 
(dotted line in Figs, to |^d). 

9 Conclusions 

Di-jet event rates have been measured in deep-inelastic scattering at small XBj (10""^ < XBj 
< 10-2) and moderate {5 < < 100 GeV^). Three different scenarios of cuts on the 
transverse momenta of the jets have been investigated: the basic symmetric requirement 
{Pt rain — 5 GeV for both jets), and additionally the asymmetric (p^min > 5 and 7GeV) 
and sum pj' (> 13 GeV) requirements. The analysis was performed in the hadronic centre 
of mass frame. The data have been corrected for detector and QED radiative effects. 
This analysis probes a region in jet phase space of small xbj and jet transverse momenta 
squared of similar size or larger than the photon virtuality {Pt/Q"^ ^1)- 

Assuming a direct or point-like photon, leading order matrix element calculations in 
combination with parton showers as an approximation of higher order effects fail com- 
pletely to describe the data. 

Adding to the leading order model additional contributions from resolving the partons 
inside the virtual photon (RAPGAP) appear to give an effective description of higher order 
effects leading to good agreement with the data for all three scenarios. A similarly good 
agreement with the di-jet event rates is observed for the colour dipole model (ARIADNE) 
with its features of gluon emission. 

Next-to-leading order calculations in assuming a point-like virtual photon provide 
a good description of the data for the scenarios with the asymmetric and the sum cut, 
except for the lowest bin in and xbj- This is improved by a NLO calculation which 
also considers contributions from resolving virtual photon structure. 
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Appendix 





symmetric (5/5 GeV) 


asymmetric (5/7 GeV) 


sum (13 GeV) 


g2 [GeV2] 


R2 


Cstat 


f^syst 


R2 


Cstat 


^syst 


R2 


Cstat 


^syst 


5-11 


0.040 


±0.002 


+0.005 
-0.004 


0.024 


±0.001 


+0.004 
-0.003 


0.022 


±0.001 


+0.003 
-0.002 


11 - 15 


0.038 


±0.002 


+0.003 
-0.002 


0.021 


±0.001 


+0.001 
-0.001 


0.021 


±0.001 


+0.002 
-0.001 


15 - 20 


0.046 


±0.002 


+0.00;-! 

-0.002 


0.028 


±0.002 


+0.003 
-0.001 


0.028 


±0.002 


+0.0()'i 
-0.002 


20 - 30 


0.050 


±0.002 


+0.006 
-0.006 


0.028 


±0.001 


+0.002 
-0.002 


0.030 


±0.002 


+0.001 
-0.002 


30 - 50 


0.067 


±0.003 


+0.006 
-0.004 


0.043 


±0.002 


+0.003 
-0.002 


0.042 


±0.002 


+0.003 
-0.005 


50 - 100 


0.093 


±0.004 


+0.006 
-0.010 


0.063 


±0.003 


+0.003 
-0.004 


0.063 


±0.003 


+0.003 
-0.003 





symmetric (5/5 GeV) 


asymmetric (5/7 GeV) 


sum (13 GcV) 


XBj 


R2 


Cstat 


^syst 


R2 


Cstat 


^syst 


R2 


0"stat 


f^syst 


10"^ - 2.5 • 10"^ 


0.046 


±0.003 


+0.006 
-0.004 


0.031 


±0.002 


-0.002 


0.029 


±0.002 


+0.001 
-0.003 


2.5 • 10"^ - 5.0 • 10"^ 


0.046 


±0.002 


+0.004 
-0.005 


0.027 


±0.002 


+0.002 
-0.001 


0.027 


±0.002 


+0.002 
-0.001 


5.0-10"''- 10-3 


0.047 


±0.002 


+0.003 
-0.005 


0.028 


±0.001 


+0.002 
-0.003 


0.027 


±0.001 


+0.002 
-0.002 


10-3-2.5- 10-3 


0.049 


±0.002 


+0.00S 
-0.004 


0.030 


±0.001 


+i).i.m 

-0.002 


0.030 


±0.001 


+0.002 
-0.003 


2.5-10-3-5.0- 10-3 


0.053 


±0.002 


+0.003 
-0.004 


0.033 


±0.002 


+IJ.i)0'2 
-0.002 


0.033 


±0.002 


+0.002 
-0.002 


5.0 - 10-3 - 10-2 


0.069 


±0.005 


+0.016 
-0.011 


0.044 


±0.003 


+0.015 
-0.007 


0.045 


±0.004 


+0.010 
-0.006 



Table 1: Di-jet rate in bins of and xbj, and statistical and systematic errors for the 
symmetric, the asymmetric, and the sum cut scenario on the pi of the two jets. The overall 
systematic error of lg°J° for the symmetric scenario and lg^J° for the asymmetric and sum 
scenario, arising from the uncertainty of the hadronic energy scale of the calorimeter and 
the uncertainty of the radiative QED corrections, is not included. 
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